Transmission and diffraction properties of a narrow slit in ideal metal 
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By solving Maxwell equations with the ideal-metal boundary conditions in the TM case, we 
have fully described the transmission and diffraction properties of a single slit regardless of its 
width. Efficiencies of the main transformation processes - transmission, diffraction, and reflection - 
are analyzed in the sub-to-few-wavelength range showing a number of sharp fundamental features. 
Close links with the case of real metal are considered. 
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Solutions to electromagnetic diffraction problems 
based on the ideal-metal boundary conditions date back 
to the works of Sommerfeld, Rayleigh, and Bethe [l|43|, 
see also [J, l5[. Such solutions are few in number, compli- 
cated, but indispensable in the subwavelength/nano op- 
tics where the usual Huygens-principle-based theory and 
intuitive approaches fail. Importantly, they link together 
the far- and near-field properties of the electromagnetic 
field including the corner singularities [a] . 

The general upsurge of interest in nanooptics gave rise 
to many research areas involving nanostructured met- 
als, such as near-field microscopy [7], metal-clad cavi- 
ties |8[ , bio-sensing y , and extraordinary light transmis- 
sion (ELT) through nano-holes [lOl lll| . Already the for- 
mulation of the corresponding problems strongly differs 
from that typical of classical optics: Instead of diffraction 
from plain obstacles f4| , one deals with funneling of light 
into (out of) apertures, with coupling of the opposite 
metal interfaces, with a local enhancement of the elec- 
tromagnetic fields. Involvement of the surface plasmons 
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in real metals further enriches the physics 

Theoretical basis of nanooptics of metals is a big issue. 
On the one hand, direct numerical methods typically map 
a tiny part of the actual space of variable parameters of 
the system with no real insight into the physics of multi- 
scale phenomena. On the other hand, they are capable of 
a dramatic enhancement of the analytical tools provid- 
ing virtually exact solutions to key physical problems. 
Reduction of complex problems to the basic elementary 
ones is a strong line of the studies, as it is known, e.g., for 
the ELT case 13l4l5|: The transmittance of a perforated 
opaque metal film can often be accurately expressed by 
the efficiencies of the elementary single- interface transfor- 
mation processes, while the film thickness trivially affects 
the positions of the Fabry-Perot resonances. 

In this letter, we present a full-scale solution to the 
single-slit problem within the paradigm of ideal metal. 
This problem is among the most basic ones in nanooptics 
of metals. Our solution reveals a wealth of subwavelength 
and near-subwavelength features which have never been 
known. While our "ideal" problem is easier than that 
for real metals (owing, e.g., to the absence of the surface 
plasmons), it is more complicated than the Sommerfeld 



problem of diffraction from a single metal wedge [2] . Fur- 
thermore, it is applicable to almost real metals. 

An important property of the slit geometry is the sur- 
vival of a single fundamental propagating mode in the 
subwavelength case for ideal and real metals. This prop- 
erty is inherent in any multiply connected aperture cross- 
section [5i[l6|, such as, e.g., a circular slit. Our results 
on the transmission properties are thus generic for a wide 
class of geometries. 
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FIG. 1. Two basic interface problems, "Out" and "In", cor- 
responding to the incidence from outside (a) and inside (b). 

Two distinct cases, "Out" and "In", for a slit of the 
half-width I are depicted in Figs, la and lb. In the case 
" Out" , a plane wave of a unit amplitude at the wave- 
length A, incident at the angle 9q, excites propagating 
modes in the slit, diffracted waves in air, and a reflected 
wave. In the case "In", a unit-amplitude propagating 
mode traveling to the interface z = reflects back and 
excites diffracted waves in air. 

In the TM case, the magnetic field amplitude has only 
a y-component H = H{x,z), and the nonzero compo- 
nents of the electric field arc Ex = (—i/ko) dH/dz and 
Ez = (i/ko) dH/dx, where ko = 27r/A. At the air-metal 
boundaries, the tangential component of E and the nor- 
mal component of H turn to zero. 

As a full set of eigenfunctions for z > we choose [5|] 
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The corresponding eigenvalues are Pi, — y/k^ — [■kv/21Y. 
The real and imaginary (3^ refer to the propagating and 
evanescent modes. For the fundamental propagating 
mode, v — Q, we have /3o — ko and ho{x) — 1. For 
/ < A/4 all other modes are evanescent. 



For the case "Out" and z ^ 0, the field H{x,z) is 
generahy presented by the expansions 

ff>=^ c,6, /i,(x)e*^-^ (1) 

1/ 



where kQj. ~ kQsmOQ, feg^ — fcocos^oi ^k = (^o ^ A:^)^'^, 
Co = 1, and c,y = 2 ioi v y^ 0. For z > and 
\x\ < ^, we have satisfied Maxweh equations and the 
boundary conditions at \x\ — I. Furthermore, we have 
bu = {H^{x,0) h^), where (..) indicates averaging over 
the sHt. In the air region, 2; < 0, we have satisfied so far 
only Maxwell equations. Here the waves with k < ko are 
propagating, while for k > ko they are evanescent. 

The amplitudes bi, and Ofc can be found if we sat- 
isfy the remaining boundary conditions: E^{x,0) — 
E^{x,0) Q{l—\x\), where Q{x) is the Heaviside step func- 
tion, and H<{x,0) = H>{x,0) for |x| < I. They are 
equivalent to the relations 
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bi, = fu, -ko^ + 9 / "fe fi^^ -k dk , 
^ J -00 

where f^^k — smc{kl + 7rz//2) -|- sinc(fcZ — tiv/2) for v = 
0, 2, . . ., fy^k = * shic{kl + 7^/2) — i sinc(fc/ — -kv (2) for 
V = 1,3,..., and sinc(..) = sin(. .)/(..). For even/odd 
v-i ]v,k is real/imaginary and even/odd in k. Combining 
Eqs. dU, we come to the set of coupled-mode equations 
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with the coupling coefficients 
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Obviously, T^^i = for the modes of different parity, i.e., 
the set ([3]) splits into two sets ~ for the even and odd 
modes. Calculating b^ from Eq. (|31) and using Eq. ^ for 
Uk, we solve completely the problem "Out". At ^o — Oi 
the amplitudes b^, are nonzero only for the even modes 
and the driving force is /^.o = 2(5^o- 

Using Eqs. ([T][2]) and the integral representations [18 1 



r(i) 



of the Hankel function ffg , one can prove lastly the 
general relation for the diffracted component of iJ<: 

H^ = -^fj^ (^'' 0) ^0'^ [ko^{x-x'Y+z^) dx'. (5) 

It generalizes the Huygens principle. The Hq function 
describes a point irradiation source placed at z' = 0, and 
the field E^{x',0), which is nonzero for |a;'| < I, serves 



as a density of the oscillating magnetic moment. The ef- 
fective irradiation source replaces indeed the real sources 
- the surface currents at the air-metal boundaries. 

The case " In" can be treated similarly. In the subwave- 
length case, which is of prime interest, the propagating 
mode incident from the inside is fundamental. Repeating 
the calculations, we come to the symmetry relations for 
the amplitudes of the excited waves/modes: a™ — —ak, 
60" = I — bo, and 6™ — —b^, ior v ^ 0, where the am- 
plitudes Ofc and bi, are taken for the normal incidence. 
Thus, the case "In" is reducible to the case "Out". 

Computation of the coefficients T^j/' and determination 
of the amplitudes b,y via truncation of Eq. ([3]) present no 
special problems, see below. However, the fundamental 
limit I -^ can be treated analytically. We have here 
Too -5- 0, f„.-k,y^ -> 2(5^0, and, correspondingly, bo = 
{H>{x, 0)) -> 2 and 6^ ^ for z/ 7^ for any 6*0. 

The dependence of the amplitude b„ on the angle 6*0 
and the ratio r = Al/X exhibits important features. The 
cases of small and large angle of incidence are different, 
see Fig. 2. For 9o < 10°, the situation is close to that for 
^0 = 0: The fundamental 0-mode dominates everywhere, 
the function l^oK^) drops from 2 to ~ 1 and then remains 
almost constant, while |6i.2,...|(^) ^ 1- For 9o > 30°, the 
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FIG. 2. The amplitudes |6o, 1,2,3! versus Alj\ for 
and 60° (b); the truncation number (--max = 20. 
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0-mode dominates only for r < 0.5. Otherwise, there is 
a strong mode competition. The opening of new propa- 
gating modes at r = 1,2,... is linked to sharp peaks of 
1^1,2,... I (j^); the odd and even peaks are comparable with 
each other tending to 1 for 60 -^ 90°. Mutual influence 
only of the modes of the same parity is evident. Only a 
few nearest modes dominate for each particular value of 
r, i.e. a selective mode excitation takes place. 

Consider now the interface characteristics. The trans- 
mission properties can be described by the efficiency 



77t(r,6'o) = ^ Cu\bu\^^efi^ I ko cos Oo 
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which is the ratio of the energy fiux through the slit to the 
flux incident onto the slit. The transmission efficiency is 
the sum of the partial contributions from all propagat- 
ing modes. For r < 1 we have rit{r,6o) = 1 60 P/ cos 6*0; 
this gives rit{0,Oo) = 4/cos0o for r ^^ 0. Similarly, we 
introduce the differential diffraction efficiency 



m{r,0,9o) = TTkocos^O \ak\^ / lcos9o 



(7) 



where 9 = arcsin(A:/A:o) is the diffraction angle. Owing to 
the symmetry properties, we have ri™{r,0) = rid{r,9,0). 
Integrating "qdir, 9, 9o) over 9, we obtain the total diffrac- 
tion efficiency r]^{r,9o). The products 2/cos0o%,d give 
the cross-sections of the corresponding transformation 
processes, and the difference Rin{r) = 1 — ?7f (?', 0) gives 
the internal- reflection coefficient for the problem "In" . In 
the limit r -J> we have r/f (0, 6*0) = and -Ri„(0) = 1. 

Owing to the non-dissipative character of the trans- 
formation processes, the amplitudes b^ obey the energy 
conservation law, which is similar to the "optical theo- 
rem" of the scattering theory. It looks especially simple 
for the normal incidence: rjt + rj^ = 2Re6o, which allows 
to express the phase of bo by r/t and 77 J for r < 1. 

Concerning the truncation procedure, it is sufficient to 
take into account 2 — 3 evanescent modes in addition to 
the propagation modes in order to achieve a sub-percent 
accuracy for ijt^d and to fulfill nicely the energy conser- 
vation law. This means that fmax ~ [4^-^] + 3. 

Fig. 3 shows the behavior of rjt and 77J when changing 
r — 4// A and 6*0. It should be considered in co-junction 
with Fig. 2. With 774(0) = 4/cos0o, the function r]t{r) is 
decreasing up to r = 1. Further increase of r results in 
bursts at the openings of new propagating modes. For 
6*0 = 0, small bursts occur only at r ~ 2, 4, . . . With in- 
creasing ^0; they appear also at r ~ 1, 3, . . . and become 
all highly pronounced and asymmetric. The asymmetry 
is due to a sharp square-root-law growth of the propagat- 
ing constants Pij{r) for r > v. The function r]^{r) grows 
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FIG. 3. The transmission efficiency (a) and the total diffrac- 
tion efficiency (b) versus 4Z/A for 60 = 0, 60°, and 80°. 

initially with a ^o-dependent slope. The further scenario 
depends on the angle 6*0. For 0o ^ 1, one sees only small 
bursts. Correspondingly, with increasing r, the reflection 
coefficient for the problem "In", Rinif) = 1 — yyf (r, 0), 
decreases from 1 to almost 0. For large 9q^ the first max- 
imum of 77 J (r, 9q\ situated deeply in the subwavelength 
range, considerably exceeds 1, and the subsequent os- 
cillations are highly pronounced and almost symmetric. 
Sharp minima of rQ(r) occur at zeros of I3ij{r), i.e., at 
the sharp maxima of \b^\{r) in Fig. 2b. 

Fig. 4 shows what happens to the angular dependence 
r]d{9) when increasing r and ^o- For ^o = Oj the quasi- 
isotropic distribution r]d{9) ~ const, which occurs for 
r < 1, transforms gradually into a central peak growing 
and narrowing for r ^ 1. After the opening of the second 
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FIG. 4. The differential diffraction efficiency for 60 = (a) 
and 45° (b). The curves refer to r = O.f, 0.5, 1.5, and 3. 



propagating mode (r > 1) this peak acquires oscillating 
tails. For ^0 > 0, the main difference in the behavior is in 
a progressive shift of the peak to the right with increasing 
r. For r 3> 1 the value of 0pcak approaches ^o and the 
diffracted waves concentrate around the reflected one. 

As seen from Fig. 4, the differential diffraction effi- 
ciency remains non-zero for 9 -^ ±90°, i.e. for ultimately 
large diffraction angles. This "grazing diffraction" grows 
remarkably and becomes strongly asymmetric with in- 
creasing ^o: 77d(90°) 3> 77d(— 90°). Depending on r, the 
maxima of 77d(±90°) correlate to those of rj^ir). The 
"grazing" diffraction is closely related to the excitation 
of the surface plasmons in real metals, see also below. 

The interface characteristics considered are sufficient 
to describe the transmission/diffraction properties of a 
subwavelength slit in an opaque film of a thickness 
d ISl [15 1 . Internal refiections of the fundamental prop- 
agating mode from the opposite interfaces lead to a se- 
quence of the Fabry-Pcrot resonances; their positions are 
controlled by the fcorf-product. At a resonance, the total- 
transmission efficiency is 774/(1 — Rin) and the intensity 
amplification factor inside the slit is rjt/ [1 — Rin)" ■ These 
quantities grow rapidly with decreasing r. At 9q ~ 0, the 
resonant transmission cross-section 21 rjt /{I — Rin) tends 
to A/2 for r ^0. 

In addition to rit^d, the spatial profiles H^(x,0) and 



E^^{x,0) show intimate features of the subwavelength 
behavior. Fig. 5 shows these profiles for the normal inci- 
dence, a few representative values of kgl, and fmax = 200. 
So high values of the truncation number are needed solely 
to resolve the corner singularities at \x\ = I. 

The profile i?<(a;, 0), see Figs. 5a and 5b, consists of 
two distinct sections (inside and outside the slit) sepa- 
rated by the inflection points at \x\ = I. In the subwave- 
length range, the function H^{x, 0) is structureless inside 
the slit and tending to the limiting value H{oo, 0) = 2 
for kol ~> 0, which corresponds to the slitless case. The 
same limiting value takes place for 9o ^ 0. This is why we 
have 60 -^ 2 and r/t — > 4/ cos^o for k^l — > 0. For kol ^ 1, 
we have already H^{x,0) ~ 1, which corresponds to a 
weakly perturbed entering into the slit and to rjt ~ 1, see 
Fig. 3a. For kol = 5, the corner value H'^{1, 0) is close to 
4/3, which corresponds to the Sommerfeld solution for 
a single metal wedge [2]. Larger values of H'^{1,0) for 



kol ^ 1 are due to the mutual influence of the wedges. 

Fig. 5c shows the near-field behavior of \E^{x, 0)| and 
\E<{x,0)\. The first function, serving as the radiation 
source in Eq. ([5]), is even in x and zero for |a:;| > I. The 
second function, which is odd in x, is not restricted to 
the slit area. At the infiection points, they both tend 
to infinity as \x — l\~^^^, which corresponds to the 90° 
corner singularities of the ideal metal [5|, ll6| . The larger 
i^max, the clearer is the singular behavior. 
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FIG. 5. Dependences 7if^(a;,0) and E^^^ixjO) in the near and 
far fields. Curves 1, 2, and 3 in a), b) are plotted for kol = 0.1, 
1, and 5. Curves 1 and 2 in c) correspond to kol — 0.5. The 
far-field profiles 1 and 2 in d) refer to kol — 0.5 and 5. 

Fig. 5d illustrates the far-field behavior. Outside the 
slit, where the characteristic spatial scale is A, we see a 
quickly establishing radiation pattern superimposed on 
the background value H^{oo,Q) = 2. The intensity 
\H'^{x, 0) — 2p decreases as l/\x\ without oscillations for 
X — I > X, as expected for the 2D case. Furthermore, it 
grows and then oscillates with increasing kol. 

Several issues are worthy of further discussion. 

- Our results show that the growth of the transmission 
efficiency r]t{r,9o) for r = Al/X — > and 9q -^ 7r/2 is 
linked to the general feature of the near-field behavior: 
The average over the slit {H{x, 0)) tends to 2 (the slitless 
limit) for r — > leading to the amplitude of the propa- 
gating mode bo = 2 for any 6*0. 

- The presence of noticeable grazing diffraction, see 
Fig. 5d, elucidates the mechanism of surface-plasmon 
generation in real metals. The radiation mechanism, per- 
taining near the slit, is almost the same for the ideal 
metal and for real metals with large negative values of 
the optical permittivity £„. The presence of the local- 
ized surface mode in the latter allows to catch the waves 
diffracted at large angles. The surface-plasmon excita- 
tion efficiency can be evaluated as r/sp — / ?7(i(7r/2), where 
/ = 2/ ^\£m\ is the numerical aperture. This estimate is 
in good agreement with [l7|. 

- To a big extent, the ideal-metal model is applicable to 
real metals with \£m\ ^ 1- The single-slit characteristics 
experience only minor changes for kol > \em\~^'^, when 



the slit width exceeds the skin depth. In many cases 
\£m\ = 10"^ — 10^, and this inequality is not restrictive. 
For l/X -^ 0, the single-slit characteristics, especially rjt, 
can experience substantial changes. 

- The corner singularities are clearly seen in the near- 
field for t'niax j^ 10^. However, they are uncoupled from 
the diffraction/transmission properties. The latter are 
linked to several lowest modes, i.e., to rough features of 
the near-field behavior. A small-radius edge rounding is 
thus not expected to produce a strong effect on rjt^d- 

In conclusion, a full-scale physical picture of the trans- 
mission, diffraction, and near-field properties of a sin- 
gle slit in ideal metal is presented. In the sub-to- 
near-subwavelength range, the transformation efficiencies 
show sharp dependences on l/X and ^o which are closely 
linked to the near-field behavior. For skin-thick films, the 
sharp subwavelength behavior leads to a strong Fabry- 
Perot enhancement of the total transmittance, and of 
light inside the slit. Oblique incidence strongly facili- 
tates the excitation of the non-zero eigenmodes leading 
to mode competition. Only few selectively excited modes 
are strongly involved in the near-subwavelength trans- 
mission/diffraction phenomena. The results obtained 
serve as a reference point in nanooptics of metals. 
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